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Abstract 

We consider pairing in a dilute system of Fermions with a short-range inter- 
action. While the theory is ill-defined for a contact interaction, the BCS equa- 
tions can be solved in the leading order of low-energy effective field theory. 
The integrals are evaluated with the dimensional regularization technique, 
giving analytic formulas relating the pairing gap, the density, and the energy 

density to the two-particle scattering length. 
PACS numbers:21.30.-x, 21.65.+f,26.60.+c 
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In the theory of fermionic matter, the expansion about the low-density limit has been 
invaluable for understanding the structure of the theory and the role of the interaction. At 
low densities, the interaction needs only be characterized by its scattering length to get 
expansions for the energy density, excitation spectrum, etc. [|TJ. However, to our knowledge 
the pairing singularity has never been incorporated into this framework. We have for example 
only the qualitative statement in ref. [JXJ] that the pairing singularity is logarithmic and 
unimportant for integrated quantities. A more quantitative statement is needed to have 
complete understanding of low-density fermionic matter. 

Another motivation for our study is the general reexamination of nuclear physics with 
effective field theory which is now taking place [0-|9j . In the effective field theory approach, 
the interaction is systematically expanded in a power series in momentum with the object 
of getting relationships between observables such that the details of the short- distance in- 
teraction need not be parameterized. We shall show here that the BCS theory of pairing 
is amenable to this approach, and the low-energy theory gives finite and analytic results. 
Within effective field theory many results can be obtained analytically opposed to the nu- 
merical treatment of potential models. In this sense our approach complements the large 
body of literature of pairing in nuclear and neutron matter that is based on potential models 

EMU- 

We consider a Fermi gas with two-fold degeneracy interacting with a short-range attrac- 
tive interaction. Examples are neutron matter or gaseous 3 He. The Hamiltonian is idealized 
to be of the form 

TT Jr f d 3 k ( | t \ Tr2 f d 3 k r d 3 k' t + 

J (2<)3 efc Kt a M + a -kd a -k,l) +gV J T^Tg J J^A^ a -Ki a -k'^k'A (1) 

where = k 2 /2m is the kinetic energy and V the volume. In effective field theory the 
contact interaction is the leading term in a derivative expansion of the many-body system. 
This limits the validity of the Hamiltonian (|I|) to the regime of long wave lengths or small 
densities. However, corrections can systematically be implemented. We have only retained 
terms in the contact interaction that are needed in the wave function. The BCS wave 
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function has the form |\&) = Hk(Uk + Vka\+a_ k i)|0); the energy is minimized with respect 
to Uk, Vf. to get the BCS equations |17[]- The equation for the pairing gap A is 

qV f d 3 k , , 

I (2) 



2 W J ^ - A) 2 + A 2 
where A is the chemical potential. The density is given in terms of these parameters by 
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Finally, the energy density of the paired state is given by 
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(4) 



(e fc -A) 2 + A 2 2 ,y (efe _ A) 2 + A 2j 

Note that the last two integrals are finite, although each integrand is a sum of terms that 
are individually divergent. 

The problem with Eq. (|j) as derived is that the contact interaction is singular in three 
dimensions. One often introduces a cutoff to make the integrals converge. However, in 
effective field theory cutoffs are not explicitly introduced. Rather, the computed observables 
are expressed directly in terms of other physical quantities. To leading order in a low-energy 
expansion of the interaction, the physical quantity is the scattering length. With the same 
Hamiltonian, the scattering length a is given by a similar divergent integral, 



mgV gV r d 3 k 1 

Ana + " 2~ J (2vr) 3 e^ 



Let us now subtract equations (0) and (|5|) to obtain 



(5) 
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2(2tt) 



d 3 k 



(6) 



(e k - A) 2 + A 2 e fc 

Notice that the integral is now convergent and so any cutoff can be taken to infinity. Fur- 
thermore, the strength of the contact interaction g, which is also an unphysical quantity, 
can be divided out. It is convenient to evaluate both terms of the integral (f]) separately 
by dimensional regularization (DR) |18|]. In DR, integrals of powers are zero so the second 
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term in the integrand in (|6]) can be dropped. The first term can be evaluated using [19 
(3.252.11), 



o 



dz m = (1 + x 2 ) Q / 2 P a (-1/VT+x^) , (7) 



where P a denotes the Legendre function. 
We write the final result in the form 



1 (i + x^Px/J-i/VTT^) (8) 



kpa 



where hp = V2mA is the Fermi momentum and x = A/A. This is our main result. Eq. (||) 
is graphed in Fig. |I[ For small values of kpa the gap is exponentially small as in the usual 
BCS theory, 

This comes about by the behavior of Pi/2(z), which has a logarithmic singularity at z — — 1 



|. Eq. (ID agrees with the result derived in ref. fT3}| . For large values of kpa, the gap is 
proportional to A, approaching A m 1.16A. 

For neutron matter the solution of eq. (^j) agrees with numerical results from potential 
models only for small values of the Fermi momentum. The large value of the scattering 
length (a = — 18.8fm) clearly limits the domain of validity of the Hamiltonian (|1|). In the 
appendix we consider pairing in the effective range approximation. This improves on the 
precision of the calculation in the low-density regime but does not enlarge the domain of 
validity. 

We complete this discussion by computing the energy density (|J) and the density @. 
The finite integrals involved are very similar to the previous one and can be evaluated using 
the same DR integral, eq. (0). The density of the BCS state is given by 



V~ A7r [+X) 



Pi /2 (-i/v / TT^ 2 ) + vTT^p 3/2 (-i/VTTx^)} , (10) 



and the energy density by 



V = - 2^^ A(1 + x2)V4 1 (1 + x2/6) Pl/2 (-VvT+^) + v/TT^ P 3/2 (-i/vT+^ 

(ii) 

For fixed density and scattering length eqs. (§) and ([10]) can be solved for the pairing gap 
and the Fermi energy. Put into eq. QTTD this yields the energy of the interacting system at 
fixed density. Fig. |2| shows a comparison with the energy of noninteracting neutrons. For 
\kpa\ ~ 1 (i.e. N/V ~ 5 x 1CT 6 fm -3 ) pairing lowers the energy about 3% confirming the 
qualitative statement that the effects of pairing on the binding are mild. 

Discussion — We now discuss the domain of validity of this low-density expansion. As 
pointed out above, the applicability of Hamiltonian (P is limited to the regime of long wave 
length \kpa\ <C 1 or small densities. The description of neutron or nuclear matter at nuclear 
densities requires the inclusion of the effective range and pions. Comparing with more 
microscopic calculations involving phenomenological potentials it appears that deviations 
from the low-density behavior are set by the scattering length. Similar considerations can 
be made for 3 He where the scattering length of the Aziz potential |2T| is large on an atomic 



scale. Many-body correlation effects will become important when kpa ml. It might be 
possible to treat them by modifying the strength of the pairing and the density of states in 
eq. ([!]). The sign would be to increase the pairing, but we have not attempted to calculate 
these effects. 

Another consideration is whether the low-density phase exists for fermionic systems 
with attractive scattering lengths. In the case of 3 He, a low-density phase could only 
be metastable at zero temperature, because there is a finite binding of the liquid phase. 
However, the metastability could be quite significant, because the minimum size for a bound 
drop is thought to be of the order of fifty particles. Another indication of the metastability 
of a low-density phase is the sound velocity in the scattering length expansion. Taking the 
first three terms, the sound velocity is positive at all densities, and thus small deviations 
from uniformity are energetically unfavorable. In the case of neutron matter, it is thought 
that pressure is always positive as a function of density, so the low-density state would be 



stable. 

In summary, we have considered the pairing in low-density Fermi systems within effective 
field theory. This model independent approach yields analytical expressions which relate the 
pairing gap, the density and the ground state energy to the scattering length. The analytical 
derivation of these results is quite interesting. 



APPENDIX 

To include the effective range we add the effective range potential 

92y2 1 / ^f^ k ~ ^HX-m^'^'.t ( 12 ) 

to the Hamiltonian ([!]). The gap equation then becomes 

V f d 3 k g + g 2 (p-k) 2 

~2J J(e k - A) 2 + 5 <13) 



and is explicitly momentum dependent. We make the quadratic ansatz A p = A + p 5 
for the momentum dependence and obtain two coupled equations that express A and 5 in 
terms of (divergent) integrals. To deal with the divergencies we observe that the integrals' 
dependence on the Fermi momentum is given by 

f d3k k2n kf +1 m . , 

/ 7o~Ts = o 2\ J n+l/2{X,y), (14) 

J (2tt) 3 . /77 _ \72 7 (a 7 u2x\2 2ir 2 \ 1 



( 27r ) 3 ^( efc - a) 2 + (a + my 

where x = A/A, y = 5k 2 F /X and J n+ i/ 2 (x, y) is the dimensionless function 



JJx,y)= / dt — (15) 
{ ^t-l) 2 + ( X + yt) 2 

In effective field theory an expansion in momenta is quite useful In what follows we 

truncate each of the gap equations to its leading order in the Fermi momentum and obtain 

Vgmk F 

— J y^y^ 

Vg 2 Amk F 

6 = 4^ Ji/2(x,y). (16) 



Obviously we have 5/A = g 2 j ' g. To make contact with low energy scattering data we expand 
the scattering amplitude 



A{p) = Vg 



(17) 



l + VgI(p) + ^p 2 + (VgI(p)) 
9 

up to quadratic order in momenta. The loop integral is 

I(p)s l[*!L ! = ?*7* tm (is) 

yF> 2 J 27rV e n - e k + in 4vr 2 J 1 - 1 + in y ' 



At low energies the scattering amplitude is given in terms of the scattering length a and the 
effective range ro 

A(p) = — [l - iap + (ar /2 - a 2 )p 2 } . (19) 
m 1 J 

Note that the divergence of the integral Ji/2(x,y) appearing in the gap equations ([16|) is 
similar to that of the loop integral I(p) appearing in the expression (|17D for the scattering 
amplitude. Thus, both divergencies may be taken care off by a a renormalization of the 
coupling constants g and g 2 - We use dimensional regularization to compute the divergent 
integrals. One obtains I(p) = —i{m / Att)p and a comparison of ( |TTD and ([19]) yields g 2 /g = 
ar /2 and g = 4na/m. Finally we have 

Ux,y) = ~-^(l + y 2 r 1/2 (\^-X' 2 P "(-*)> ( 2 °) 

sin 7ia y 1 + y z I 



where z — (1 — xy)/ J (1 + x 2 )(l + y 2 ). This yields the final results 



' /- ■ 2\-l/2 / 1 + x 



2\ 1/4 



y = °Yk 2 F x. (21) 

Note that these equations add corrections of the order ~ k 2 F ar Q to the gap equation 
These corrections are small only in the low-density regime kpa 1. For a description of 
neutron matter (a = — 18.8fm ro = 2.75fm) at larger densities, at least the inclusion of 
pions seems to be necessary. Note also, that the gap equations fl2~T| ) become singular for 
k 2 F — > —2/ar (i.e. y — > —x) due to the logarithmic singularity of the Legendre function 



for z — > — 1. This behavior results from the quadratic approximation for the interaction 
potential and the truncations in the gap equation. It is related to the change in sign of 



the truncated potential at kp = y—2/aro |TJj]. Again, the introduction of pions or higher 
potential terms seem to be necessary to alter this behavior. 
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FIGURES 
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FIG. 1. Energy gap A/A as a function of fcpo 
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FIG. 2. Energy of interacting neutron matter normalized to the energy Eq of noninter acting 
neutrons as a function of the density. 
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